A spaceship with a thruster - one body, one force 
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A spaceship with one thruster producing a constant magnitude force is analyzed for various 
initial conditions. This elementary problem, with one object acted upon by one force, has value as a 
challenge to one's physical intuition and in demonstrating the benefits and limitations of dimensional 
analysis. In addition, the problem can serve to introduce a student to special functions, provide a 
mechanical model for Fresnel integrals and the associated Cornu spiral, or be used as an example in 
a numerical methods course. The problem has some interesting and perhaps unexpected features. 



I. INTRODUCTION 

A problem involving one constant magnitude force act- 
ing on one body leads to interesting motion and is useful 
as a teaching tool when discussing physical intuition and 
dimensional analysis. The problem can be solved with 
elementary mechanics with the solution expressed in the 
form of integrals associated with known special functions. 
We first state the problem and then proceed with its solu- 
tion, though we invite the reader to ponder the problem 
before reading the solution and the remainder of the ar- 
ticle. After the solution is presented we discuss what 
one could have ascertained from an astute application of 
physical intuition. Next, dimensional analysis is applied 
to the problem. We hope instructors will find the richness 
of such a simple to state problem to be of value in teach- 
ing and in challenging students. In addition, the problem 
can be attacked analytically to a large degree and thus 
comparison of analytic (and asymptotic expressions) to 
numerical solutions could be instructive in a numerical 
methods course. 



II. THE PROBLEM 

Imagine a spaceship with one thruster positioned a dis- 
tance R from the center of mass as depicted in Fig. ^ 
Initially the ship is at rest. At time t = 0, the thruster 
is fired and produces a constant tangential force, F. De- 
scribe the motion of the ship. What path does the cen- 
ter of mass move along? Assume special relativity is not 
needed and that the mass of the spaceship/thruster com- 
bination does not change. What else can one intuit about 
the motion? The solution follows immediately so we sug- 
gest the reader formulate opinions about the motion be- 
fore proceeding. 




FIG. 1; Position of spaceship/thruster combination at time 
t = 0. 



Spaceship is now rotating 

F 



Path of the 
center of mass 




Thruster delivering 
constant magnitude force, F 



FIG. 2: Qualitative sketch of the position of the space- 
ship/thruster combination shortly after t = 0. The ship is 
rotating and translating. 



A. The solution to the problem 

Three equations, one from the torque of the thruster, 
and the other two from the x and y components of the 
force of the thruster, define the motion. Let 9 be the 
angle the thruster has rotated about its center of mass 
since time t = 0. Initially, the ship will move in the x 



direction and then upward into the first quadrant as it 
begins to rotate as shown in Fig. [3 

From the torque on the system about the center of 
mass, the rotational analog of Newton's second law (r — 
16) requires 

FR^cmR^e (1) 

where the moment of inertia is / = c mR^ , where c is a 



dimensionless constant that depends on the distribution 
of mass (e.g. i for a disk, but in general any positive 
number), m is the mass of the ship, R is the distance 
from the center of mass to the thruster, and 9 is the 
second derivative of 9{t) with respect to time. Solving 
for the angle as a function of time we have 
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assuming 9{Q) = 9{0) = 0. The case of ^(0) ^ will be 
explored later in the paper. 
Newton's second law requires 
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where x and y are the coordinates of the center of mass. 
Substituting for 9 and integrating these equations gives 
velocity components: 
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assuming Vx{Q) = Vy{Q) — 0. These integrals are well 
studied and are called Fresnel integralai. Their evalua- 
tion is aided by a plot called the Cornu Spiral, which is 
shown in Fig. |3| Note the analysis required to this point 
was within the level of the typical elementary calculus- 
based physics course, though the integrals have led to 
special functions. 



B. The motion of the center of mass and the 
Cornu Spiral 

Examining the Fresnel integrals for the velocity com- 
ponents reveals much about the motion of the center of 
mass. In the limit t — > cx), each component of veloc- 
ity goes to yJiTcFR/Am. Therefore, the center of mass 
moves off at a 45 degree angle with respect to the x- 
axis as time approaches infinity. The Cornu Spiral, most 
commonly associated with the problem of diffraction of 
a rectangular aperture^, represents these two integrals 
graphically. Fig. O is the Cornu Spiral as it is often 
displayedi^. A point on this spiral (when multiplied by 
the factor ^JircFR/m) represents the x and y compo- 
nents of velocity at a point in time. A line drawn from 
the origin to a point on the spiral represents the instanta- 
neous direction of the motion of the center of mass. Thus 
from this plot we see the center of mass is always in the 
first quadrant since v^ and Vy are always positive. And 
since v^ is positive at all times then the plot of motion 
in the y versus x plane is single valued. 
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FIG. 3: The Cornu Spiral as it is often shownlia. The Fresnel 
sine integral, S{z) = J sm{^t^)dt, is plotted on the y-axis 
versus the Fresnel cosine integral, C{z) — J^ cos{^t^)dt, on 
the a;-axis. The parameter z, proportional to time for our 
problem, runs along the spiral with a tick every 0.2 stopping 
at 3.0. For the spaceship problem the axes are proportional 
to the X and y velocity components. The central portion 
of the spiral near (0.5, 0.5) is not plotted for clarity. The 
three straight lines on the plot represent the velocities with 
the steepest direction, the terminal speed, and the maximum 
speed, respectively in a clockwise fashion. 



Fig. 121 also reveals a maximum speed (depicted by the 
longest straight line) which is approximately equal to 
1.3422wtx), where Voa (straight line to center of spiral) rep- 
resents the speed of the center of mass as time approaches 
infinity. The plot also shows a maximum angle with re- 
spect to the X-axis for the trajectory approximately equal 
to 60.466°, independent of any other parameters. Thus 
if the thruster is fired for an appropriate finite period 
of time, one could obtain a trajectory anywhere between 
zero and 60.466° with respect to the x-axis, or a terminal 
velocity anywhere between zero and 1.3422tioo. 



C. The path of the center of mass 

Integrating the velocity components, equations El and 
El gives the position of the center of mass: 
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(8) 



where we have assumed x(0) = y(0) = 0. These inte- 
grals can be evaluated numerically. Interestingly, using 



Y 






7 
6 
5 
4 
3 
2 
1 






^ 






■^ 















1 



3 



7 



FIG. 4: The path of the center of mass of the spaceship. For 
actual distances the axes should be multiplied by cR. The 
shape of this path is universal and doesn't depend on other 
parameters save for the scaling factor cR. 



integration by parts the position can also be expressed 
analytically in terms of equations [S] and [^ the compo- 
nents of velocity. We find 



x{t) ~ —ci? sin( 



Ft' 
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(9) 



(10) 



Thus the path of motion can be studied analytically 
through the Fresnel special functions and its associated 
Cornu spiral. 

A plot of the motion is shown in Fig. 21 The shape 
of the path is universal regardless of parameter values 
though distances are scaled by the factor c R. Note that 
the asymptotic trajectory as projected back towards the 
origin does not pass through the origin but has a non-zero 
x-intercept. Analytic analysis of the asymptotic (t — > oo) 
forms of equations El and ^1 show this intercept occurs 
X = cR. Curiously, while the thruster delivers equal 
X and y components of impulse (change in momentum) 
as time approaches infinity to the center of mass, there 
is an asymmetry in the displacement as shown by this 
intercept (due to Vx being greater than Vy initially). For 
actual spacecraft maneuvers we see a single thruster is 
not a very practical configuration. 



III. PHYSICAL INTUITION 

It has been our experience that the majority of stu- 
dents and faculty alike have difficulty intuiting the mo- 
tion of the center of mass. Physical intuition is not a 



well defined term. An interesting recent book entitled 
Seeking Ultimates: An Intuitive Guide to Physics^ states 
intuition is something for a student "to absorb in their 
bones." The dictionary^ defines intuition as 

la) the act or faculty of knowing without the 
use of rational processes; immediate cognition 
b) knowledge acquired by use of this faculty. 
2.) acute insight 

We feel a definition of "physical intuition" requires more. 
A recent article by Singh^ agrees that physical intuition 
is difficult to define but offers these words: 

Cognitive theory suggests that those with 
good intuition can effectively pattern-match 
or map a given problem onto situations with 
which they have experience. 

These words provide a suitable footing for the term be- 
cause below we relate the problem at hand to a more 
common problem, which most physicists have had expe- 
rience with during the course of their education. Perhaps 
"absorb in their bones" is on the mark if interpreted as 
absorbing a number of standard problems to provide a 
bank with which to pattern match. 



A. Center of mass has a terminal velocity 

The simplest idea is that as the object spins faster and 
faster the impulse to the center of mass over a single 
revolution must tend to zero. Therefore the change in 
linear momentum tends to zero and thus the notion of 
a terminal velocity for the center of mass is reasonable 
(though not guaranteed, the harmonic series tends to zero 
but it's sum does not). 

The first half of a revolution takes longer than the sec- 
ond and thus it must be the case that the impulse is 
always positive in the y direction for any time and the 
motion is confined to the upper half plane. In addition, 
y plotted as a function of time is monotonically increas- 
ing. One may be tempted to draw similar conclusions for 
the X direction but here things are trickier, especially for 
whether a plot of x versus time is monotonically increas- 
ing. To see this consider a slightly different problem. 



B. An alternate spaceship problem 

If the spaceship's thruster had acted through its center 
of mass and had a rotation rate given by ujo, as pictured 
in Fig.|Sl then we could say something about the x com- 
ponent of velocity. During the first quarter of rotation 
the X component of acceleration is positive and the x 
component of velocity goes from to some maximum. 
During the second quarter of rotation the x component 
of acceleration is negative and the symmetry of the ap- 
plied force dictates that this acceleration will reduce the 



Thruster acts through center 
and is fixed to the rotating 
spacecraft 




FIG. 5: Sketch of the initial position for the alternate problem 
with the thruster acting through the center of mass but with 
the ship initially rotating at angular velocity uOo- 



X component of velocity to zero. During the last half of 
the rotation the x component of velocity will be negative 
and the symmetry of the kinematics would return the x 
component of the center of mass to a; = 0. Then, as far 
as the X direction is concerned, the whole thing starts 
over again. The overall motion, assuming 9{t) = LJot, is 
a cycloid, reminiscent of the motion of a charged par- 
ticle starting at rest in orthogonal uniform electric and 
magnetic fields. 

The velocity components are 

F 

Vx{t) = 3— sin(a;ot) (11) 
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The path is depicted in Fig. and the shape of the 
path is also universal, though the axes are scaled by the 
factor F/tow^. One possible mistake is to confuse con- 
stant rotation with uniform circular motion. But uniform 
circular motion is not a correct analogy since the force 
of the thruster is not, in general, perpendicular to the 
velocity of the center of mass. 



C. Spaceship is stuck in the first quadrant 

Returning to the original problem, the rotation rate is 
not constant, but increases. As such we would expect 
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FIG. 6: The cycloid motion of the center of mass for the 
alternate problem. The shape of the path is universal. The 
axes are in units of F/miol . 



the particle to never return to x = since the time spent 
in each rotation thrusting with a positive x component 
will be longer than the time spent with a negative x com- 
ponent. Thus the spaceship is doomed to remain in the 
first quadrant for all its travels contrary to a common 
misconception that the spaceship may move in some sort 
of spiral around the origin. 

As mentioned earlier, we note the actual path of the 
center of mass as described by a function y{x) is single- 
valued, meaning physically that the x component of ve- 
locity (as well as the y velocity component) is always 
positive. However, had there been an initial rotation, 
uj^cmR/ F — 1 for example, then there would have been 
a negative x component of velocity during the first ro- 
tation and thus y{x) would have been double valued for 
some X values as shown in Fig. {7\ The situation of a 
thruster with initial rotation is discussed in detail below. 

It is hoped the above discussion sheds some light on 
why the 45-degree asymptotic path, the non-zero x- 
intercept, and the single- valued nature of y{x) are dif- 
ficult to intuit, even in hindsight. They depend on the 
value of integrals that arc not intuitive (without the aid 
of the Cornu spiral or some other such device) . 

We did succeed in providing an "intuitive" explanation 
to explain that the path of motion is all in the first quad- 
rant by comparing to an alternative known elementary 
problem. And, we intuited the notion of a terminal ve- 
locity and thus the asymptotic path for large times is a 
straight line. 

In the interest of full disclosure, we add that our ini- 
tial thoughts on the motion weren't always right, and 
we wrote this section with the benefit of hindsight from 
solving the equations of motion. 
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FIG. 7: Path of motion of the center of mass for a situation 
with initial rotation. Again, as in Fig. 2] the shape of this 
path is fixed for a given value of lolcmR/ F , which in this 
case is 1. The axes are in units of cR. Note the x component 
of velocity is negative for a period during the first rotation 
(near x = y — 1). 



IV. DIMENSIONAL ANALYSIS 



TABLE I; Quantities of interest, dimensional estimates, and 
actual values for the problem with the thruster at the edge of 
the spaceship. 



Quantity 




Estimates 


Actual 
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mass 
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st rotation) 


R 

m 


- 


time 


y/mR/F 
\jmRIF 
l/FR/m 

R 

R 


- 


^Ist rotation 
Voc 

X intercept of a 
displacement (fi 


^jA-KcmR/F 

^■KcFRjIm 

cR 

~ 3.75cR 



TABLE II: Quantities of interest, dimensional estimates, and 
actual values for the alternate problem with the thruster act- 
ing through the center of mass but with initial rotation, ujo- 



Quantity 


Estimates 


Actual 


length 


F/muja 


- 


mass 


m 


- 


time 


I/OJO 


- 


^Ist rotation 


l/i^o 


2-n/uJo 


Vdrift 


F/mUo 


F/muJo 


maximum x 


F/mu)l 


2F/mLol 


displacement (tist rotation 


) F/rauol 


2-KF/mujl 



Dimensional analysis has been discussed, for example, 
in association with models and data utilizing the simple 
pendulum as an example^, in a simple experiment in- 
volving the flow of sand^, and in the error analysis of a 
falling body'*. This problem lends itself to dimensional 
analysis, the most interesting example being the terminal 
velocity. The characteristic mass is m, length is i?, and 
time is ymR/F . A fourth parameter for the problem 
is actually dimensionless, the parameter c from the form 
factor of the moment of inertia. Though dimensionless 
one can usually predict whether a quantity should in- 
crease or decrease as a function of c, though the power 
of the dependance on c is unobtainable by such analysis. 

TableQlhst a few quantities of possible interest, such as 
terminal velocity [voo ) and the x intercept of the asymp- 
totic path, along with the actual value and a dimensional 
estimate. Note all numerical prefactors of the estimates 
are within a factor of ten of the actual prefactor. 



A. Dimensional analysis of the alternate problem 

To physically understand the motion we introduced the 
alternate problem of a spaceship initially rotating with 
thrust acting through the center of mass. This effec- 
tively eliminates the radius, i?, from the problem (since 
no torque is available the rotation rate will not change), 
but it introduced a new parameter, the initial rotation 



rate, uJo- For this problem the characteristic mass is ttt., 
length is F/mw^, and time is I/wq. 

Table im is analogous to TablelUfor this alternate prob- 
lem. Note again all numerical prefactors of the estimates 
are within a factor of ten of the actual prefactor. 



B. Original problem with initial rotation 

If the original problem had been initially rotating then 
there would have been two length scales, two time scales, 
and even two mass scales. The second mass scale would 
be given by F/Ruj"^. With two sets of characteristic 
scales, dimensional analysis is of less value because there 
are an infinite number of ways to construct quantities 
of interest. For example, let li,mi,ti be a characteris- 
tic length, mass and time respectively, and let ^21 "^2; ^2 
be a second set. Suppose we're curious about a velocity. 
Obvious possibilities are h/ti and l2/t2- But 
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t2' 
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(15) 



are examples of other possibilities. 

Consider the initial problem but now allow an initial 
rotation rate as well. The velocity components would be 
given by the integrals: 
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(17) 



Our intuition says the terminal velocity should decrease 
as uJo increases (for positive Wo, i-C. in the direction of 
the applied torque). Dimensional analysis for a velocity 
reveals ambiguities such as: 



F2 



i^3/2 



F 



(18) 



The new velocity integrals can still be interpreted with 
the aid of the Cornu spiral. By completing the square of 
the arguments of the trigonometric functions, an initial 
rotation can be shown to be a shift of ^Juj'^cmR/'KF 
along the spiral and a rotation of axis by an angle of 
u!^cmR/2F as shown in Fig.|Hl where the shifted axes are 
placed for an initial rotation satisfying uj^cmR/F — 1. 
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FIG. 9: The path of the center of mass for a situation with 
initial rotation opposite the applied torque of the thruster, 
assuming ulcmR/ F — 25. Axes are in units of the factor cR. 
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(20) 



Each term in the expansions are further examples of the 
ambiguity in constructing velocities with two character- 
istic scales. The asymptotic trajectory (angle from the 
X-axis) approaches 90 degrees, since tan(0) = (vy/vx) ~ 
Lo^cRm/F a,s LOo ^ OO. 



FIG. 8: The Cornu spiral is shown with a new axis placed 
along the spiral to account for an initial rotation satisfying 
LOoCmR/F = 1. Note the trajectory angle of the line which 
is proportional to Woo (line from the new origin to the point 
|, ^) is parallel to the asymptotic trajectory in Fig.|7|above. 
The ticks placed along the spiral represent other origins cor- 
responding to LOoCmR/F = 4,9 and 16. As in Fig. |31 actual 
velocities are obtained by multiplying the axes by the factor 
^y^TcFR/m. 



With the aid of this shifted axis, we see the termi- 
nal velocity should indeed get smaller and also the angle 
of the trajectory should increase. Also, both Vx and Vy 
should approach zero as lUq approaches infinity. However, 
since the trajectory tends to 90 degrees from the x-axis 
as Wo — > OO we note v^ and Vy cannot tend to zero with 
the same dependence on cjo- 

In fact, it can be shown that as Wo — > oo 
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PATH OF MOTION WITH A NEGATIVE 
INITIAL ROTATION 



Since we have just generalized the original problem to 
include a non-zero initial rotation aligned with the ap- 
plied torque, it is interesting to consider a negative initial 
rotation, i.e. initially spinning opposite the direction of 
the appHed torque. Fig. displays the path of the center 
of mass for the situation oj'^cmR/ F = 25 and with ujo 
being negative (i.e. opposite the direction of the torque). 
Notice that the displacement vector, for this case, sweeps 
a polar angle somewhere between 270 and 360 degrees. 
This raises questions: What is the maximum this angle 
could be? Could the spaceship spiral around the origin, 
with an appropriate initial rotation, as some incorrectly 
suggest for the original problem with no initial rotation? 
Our explorations reveal that with an appropriate choice 
of a negative Uo the asymptotic path can be any com- 
pass heading in the full 360 degree range of possibilities. 
Fig. llOlshows four such possibilities associated with four 
different initial rotation rates. 
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FIG. 10: The path of the center of mass for a situation with 
initial rotation opposite the applied torque of the thruster, 
with ujo = —3.2, —3.6, —3.8 and —4.2 respectively in units of 
^J F/cmR (left to right then top to bottom). Axes are in 
units of the factor cR. 



A. Can the center of mass circle the origin? No 

As before, using integration by parts the components 
of position can be expressed in terms of the velocity com- 
ponents. We find: 



x((t) = —ci? sin( 



2cmR 



■ UJot) + {t ■ 



LOoCmR 



)vAt) (21) 



y{t) = ci?[cos( 



Ft' 



2cmR 



+Woi)-l] + (t + ^^)«,(i) (22) 

where Vx (t) and Vy (t) are those given by equations lid 
and 1 1 71 respectively. For the spaceship to circle the origin 
a necessary, but not sufficient, condition is that Vy — 
while y > 0. Examining equation] 
then 



we see that iivy = 



y\vy=Q = cR[cos{ 



Fe 
2cmR 
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(23) 



must be positive if we are to meet this condition. But 
this cannot ever be true since the most the cosine term 
could be is 1. Therefore once y > the spaceship will 
remain in the upper half plane. Note the above proof 
does not require any properties of the Fresnel functions, 
just integration by parts. 



B. What is the displacement vector's maximum 
polar angle? 

Numerically, we have determined the maximum angle 
the displacement vector can sweep while going around 
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FIG. 11: The polar angle swept by the displacement vector 
(in degrees) from i = Otoi = ooasa function of the initial 
rotation rate uJo in units of ^jF/crnR. The discontinuities 
correspond to paths that intersect the origin. Thus, qualita- 
tively, at the discontinuity the path went from one like that 
in Fig. |5] where it was going around the origin to one like that 
in the top right of Fig. llOl where it cut back below the origin. 
The maximum angles referred to in the text are defined as 
negative angles on this plot. 



the origin is approximately 319.52° which occurs when 
the initial rotation rate is approximately uJo ~ —3.54 in 
units of y/F/cmR. At approximately ujo = —5.01, which 
is near that depicted in Fig.|51the path again approaches 
intersection with the origin with a corresponding maxi- 
mum polar angle for the displacement vector of approxi- 
mately 318.27°. There are infinitely many more of these 
pairs; the list begins like this 
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i -6.13996 
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a 317.327 



(24) 



The maximum polar angle appears to continue to de- 
crease. Thus 319.52° appears to be the approxi- 
mate maximum regarding encircling the origin. A plot 
of the polar angle swept by the displacement vector 
{Q displacement) from i = to t = oo is showii vcrsus the 
initial rotation (wo) in Fig. II II in units of ^J F / cmR. 



VI. CONCLUSION 

As one object (see appendix), one force problems go, 
this one may rival the simple harmonic oscillator for its 
richness. The problem has utility in introducing a stu- 
dent to special functions and handbooks such as Abro- 
mowitz and Stegun^ . It also provides a mechanical model 
for thinking about Fresnel integrals. It could be used in 
a numerical methods course where comparison between 
analytic (and analytic asymptotic) expressions versus nu- 
merical techniques could be performed. The dimensional 
analysis applied to this problem is useful for many other 
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problems. For example, projectile motion possesses a 
universal path shape, a parabola, characterized by the di- 
mensionless parameter the launch angle with the length 
scale set by v^j g where Vo is the initial velocity and g 
the acceleration due to gravity. When presenting a new 
problem to a student, a good question to ask is to try 
to sort out how many sets of scales does the problem en- 
compass and what can dimensional arguments say about 
the answers to any questions posed. Finally, the problem 
is a challenging test of physical intuition and it can be of 
interest to the teacher and student alike to think about 
just what is meant by such a term as "physical intuition" 
and how would one go about improving it. 



Integrating the a;-component, for example, over one 
revolution (from to I-k) should be equivalent to evalu- 
ating Eq. [SI multiplied by to from i = to ^jAmcraRjF^ 
i.e. 



F 



V' 



cos(- —)dt 



2tt 



FcmR 
20 



cos{e)de 

(A.2) 

With the substitution 6 = Ft^ /2cmr this is shown 
to be true. In fact an alternative representation of the 
Fresnel cosine integral isi: 
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APPENDIX: THE OTHER OBJECT(S) - 
MOMENTUM CONSERVATION 

The title stated there is only one body in this prob- 
lem and indeed the spaceship with its attached thruster 
has been our focus. But momentum conservation sug- 
gests this cannot be the only thing in our universe. The 
thruster must be emitting something (perhaps a photon) 
that carries momentum (and also energy and angular mo- 
mentum). The momentum is carried away in all direc- 
tions since the spaceship rotates. The magnitude of the 
instantaneous impulse imparted by the thruster, is Fdt. 
The impulse per angular bin from 6 to 6-\-d9 as & function 
of 9 using Eq. |21 is then: 



Fdt _ FcmR 
~df^ V 29 ■ 



(A.l) 



A plot of this impulse density over the first cycle is shown 
in Fig. [121 



' , 2x , 1 /■" cos(a:) 



dx. 



(A.3) 
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FIG. 12: The distribution of Fdt/d9 during the first rotation 
as a function of the rotation angle. The average over the first 
cycle occurs at '^ 50.77° (~ 0.886 radians and denoted by a 
solid grid line), which is near that of the asymptotic path of 
45 degrees. 
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